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Abstract 

The paper deals with a theoretical investigation of the peristaltic transport of a phys- 
iological fluid in a porous asymmetric channel under the action of a magnetic field. The 
stream function, pressure gradient and axial velocity are studied by using appropriate an- 
alytical and numerical techniques. Effects of different physical parameters such as perme- 
ability, phase difference, wave amplitude and magnetic parameter on the velocity, pumping 
characteristics, streamline pattern and trapping are investigated with particular emphasis. 
The computational results are presented in graphical form. The results are found to be in 
perfect agreement with those of a previous study carried out for a non-porous channel in 
the absence of a magnetic field. 
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1 Introduction 

Peristaltic transport is a natural mechanism of pumping most physiological fluids induced by a 
progressive wave of area contraction or expansion along the length of the boundary of a fluid- 
filled distensible tube. Different physiological phenomena, such as the flow of urine from kidney 
to the bladder through ureters as well as the transport of food material through the digestive 
tract, semen in the vas differences, fluids in lymphatic vessels, bile from gall bladder into the 
duodenum, movement of spermatozoa in the ductus efferentes of the male reproductive tract and 
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cervical canal and the transport of ovum in the fallopian tube, the movement of cilia, circulation 
of blood in small blood vessels, take place by the peristalsis mechanism. This phenomenon is 
also observed in the propulsion of some industrial fluids. 



Nomenclat ure 




ai,6i 


Wave amplitudes 


Bo 


Intensity of the magnetic field 


di,d2 


Width of the channel 


Hi, H2 


Vertical displacements of the upper and lower walls respectively 


M 


Magnetic parameter 


P 


Fluid pressure 


Q 


Flux at any axial station 


R 


Reynolds number of the fluid 


t 


Time 


u,v 


Velocity components in X,Y directions respectively 


X,Y 


Rectangular Cartesian co-ordinates 


a 


Function defined in equation (12) 


S 


Wave number 


Ap 


Pressure difference between the channel ends 


X 


Wave length of the traveUing wave motion of the wall 


1^ 


Viscosity of the fluid 


V 


Kinematic viscosity of the fluid 





Phase difference 


P 


Density of the fluid 




Stream function 


a 


Electrical conductivity of the fluid 



This mechanism flnds applications in roller and flnger pumps, heart-lung machine, blood 
pump machine, dialysis machine and also noxious fluid transport in nuclear industries. Some 
of the recent important theoretical studies on peristaltic transport include those by Usha and 
Rao [1], Misra and Pandey [2, 3, 4, 5, 6, 7], Mishra and Rao [8, 9], Mekheimer [10], Aflfl [11], 
Eytan et al. [12], Eldabe et al. [13], Hakeem et al. [14] as well as Hayat and Ali [15]. These 
studies have been carried out either by using lubrication theory neglecting the fluid inertia and 
wall curvature without any limitation of wave amplitude or by considering small peristaltic wave 
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amplitude with arbitrary Reynolds number. Using the assumptions of thin shell and lubrication 
theories, Carew and Pedley [16] put forward a mathematical analysis for the study of peristaltic 
pumping in the ureter. Antanovskii and Ramkissoon [17] also took the help of lubrication theory 
in order to study the peristaltic transport of a compressible viscous fluid through a pipe when 
the pressure drop changes with time, by taking into account the wall deformation of the pipe. 
Lee and Fung [18] studied theoretically the peristaltic phenomenon observed in blood flow, while 
Lee [19, 20] reported his experimental studies on peristalsis observed in large intestines. 
Taylor [21] conducted a theoretical study on asymmetric wave propagation in waving sheets with 
an object to put forward an explanation to the mechanical interaction between neighbouring 
spermatozoa. By considering a biharmonic stream function in the form of a Fourier series and 
retaining only the starting term of the stream function series, he made an attempt to simulate 
cases for symmetric /asymmetric wave propagation. The boundary integral method along with 
a suitable numerical technique was employed by Pozrikidis [22] to study the fluid dynamics 
involved in the problem of peristaltic transport in an asymmetric channel under Stokes flow 
conditions. A mathematical analysis for peristaltic transport of an incompressible Newtonian 
fluid was carried out by Mishra and Rao [8] as well as Elshehaway et al.[23]. In both these 
studies, channel flow was considered where the channel is asymmetric. These two studies arc 
similar in nature, however unlike [8], the channel was considered to be porous in [23]. The work 
reported in [8] was extended by Reddy et al. [24] to the propagation of peristaltic transport of a 
power-law fluid in an asymmetric channel. This theoretical study was based on a similar study 
conducted earlier by Reddy et al. [25] for the propagation of symmetric peristaltic waves in a 
Newtonian fluid model. 

Recent physiological researches (Cf. Vries et al.[26]) indicate that uterine peristalsis result- 
ing from myometrial contraction can take place in both symmetric and asymmetric directions. 
In different studies pertaining to the gastrointestinal tract, intra-pleural membranes, capillary 
walls, human lung, bile duct, gall bladder with stones and small blood vessels, flow in porous 
tubes and deformable porous layers have been examined by various authors. It is known that 
the primary function of gastrointestinal tract consists of absorbing nutrients from the mix of 
food and liquid that passes through the tract. As mentioned by Keener and Sncyd [27], gastroin- 
testinal tract is surrounded by a number of heavily innervated muscle layers which arc smooth 
muscles consisting of many folds. There exist pores in the junctions between them, although 
the junctions are tight. Also, Bergel [28] observed that the capillary walls are surrounded by 
flattened endothelial cell layers which are porous. 



3 



Currently studies on peristaltic motion in magneto hydrodynamic (MHD) flows of electrically 
conducting physiological fluids have become a subject of growing interest for researchers. This 
is due to the fact that such studies are useful particularly for having a proper understanding of 
the functioning of different machines used by clinicians for pumping blood and magnetic reso- 
nance imaging (MRI). The related studies also contribute to updating the information that is 
required to operate machines like MHD peristaltic compressors. As pointed out by Agrawal and 
Anwaruddin [29], theoretical researches with an aim to explore the effect of a magnetic field on 
the ffow of blood in atherosclerotic vessels also find application in a blood pump used by cardiac 
surgeons during the surgical procedure. Investigations related to the study of the effect of a 
magnetic field on the movement of blood indicate that if we apply a frequency pulsing magnetic 
field of low intensity in a controlled manner, it is possible to modify the behaviour of cells and 
tissues. Li et al. [30] observed that an impulsive magnetic field can be used as a theraptic means 
to treat patients who have stone fragments in their urinary tract. 

In view of all the aforesaid observations and remarks, the present theoretical study has been 
designed in such a manner that it can explore a variety of information that can enrich the in- 
formation bank required for magneto therapy, which consists of applying a magnetic field for 
the treatment of various diseases. From the physiological fluid dynamics point of view, it is a 
study of the peristaltic transport of a physiological fluid in a porous asymmetric channel under 
the action of an external magnetic fleld. The peristaltic wave train on the channel wall has 
been considered to have different amplitudes and phases in view of the variation in the channel 
width, wave amplitude and phase difference. It is believed that the results presented here will 
flnd prospective application in the study of various fluid mechanical problems associated with 
gastrointestinal tract, intra pleural membranes, capillary walls and small blood vessels. 

2 The Problem and its Solution 

Let us consider the peristaltic motion of an incompressible viscous fluid through a porous medium 
channel under the action of a magnetic field. Let Y = Hi and F = i/2 be respectively the upper 
and lower boundaries of the channel. The medium is considered to be induced by a sinusoidal 
wave train propagating with a constant speed c along the channel wall, such that 

Y = Hi = di + aicos{^{X — ct)) on the upper boundary 
and 

Y — H2 — —d2 — biCos{^{X — ct) + 0) on the lower boundary. 
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where oi, 61 are the wave amphtudes , A is the wave length, di + 6,2 the width of the channel, 
the phase difference (0 < < tt) of the waves. = corresponds to the symmetric channel 
with waves out of phase and (f) — n represents the waves in phase; quantities ai, 61, di, d2, and 
(f) are such that 

+ 6^ + 2aibiCos{(P) < {di + d2f (1) 



The equations that govern the fluid motion in a porous medium under the influence of a magnetic 
field are: 

dU dV ^ 

ax + W = ° 

dV ,dV ,,dV IdP fe'V d'V\ i^. , 

where = ^, k is the permeability parameter, Bq the intensity of the magnetic field acting 
along the y-axis and the induced magnetic field is assumed to be negligible. 
Considering a wave frame {x,y) moving with a velocity c away from a fixed frame {X,Y), let 
us the transformations 



x^X-ct, y^Y, u^U-c, v^V, p{x) = P{X, t) (5) 

in which {u, v) and {U, V) arc the velocity components, p and P stand for pressure in wave and 
fixed frame of reference respectively. 
If ip be the stream function, we have 

Now eliminating P between (3) and (4), we obtain the governing equation expressed in terms of 
■0. Thus the equation governing the steady fiow of the fiuid in the wave frame of reference can 
be written in the form 



'4^y'4^yyx '4^x'4^yyy ~l~ '4^y'4^xxx '4^x'4^xxy ^(.^yyyy ~^ '^'4^xxyy ~l~ '4^xxxx) 



J^Ay - - -yAy (6) 
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In the sequal, we shall make use of the following non-dimensional variables: 

X _ y _ u _ V . di _ djp - ct . Hi H2 , d2 

X = T-' y=l~^ u=-: '^=-?> <^=-r' — 7' ^=T' ^1 ^ TT' "2 = -^-, d=—, 
A di c CO A iicA A d\ di di 

(^1 , bi ^ cdi T '0 r ^ , , O" ^9 ,9/ X 

« =T,' « = - ^ = ;^' = ^B.^'iSO 2) (7) 

Henceforward we shall use the non-dimensional variables defined in (7). If we now drop bars 
over the symbols, the equation governing the flow of the fluid reads 



RS ^'4^y'^yyx '4^x'^yyy ~l~ ^ {,'4^y'^xxx '^x'^xxy)^ ^yyvv ~^ ^xxyy ~l~ ^ "^xxxx 

1 5^ 

- J^i^yy - J^i^xx - Mipyy (8) 

Boundary Conditions: 

Considering the physics of the problem, the boundary conditions may be put mathematically as 



q 

(0 ip = -(ity = hi = l + a cos{2t:x) 

■0 = — ^ at y — h2 — —d — h cos{2'kx -\- 0) 

in) — = -1 at hi, /i2 (9) 
dy 

where q refers to the flux in the wave frame and a, b, d, are such that 

a^ + b'^ + 2abcos{(j)) < (1 + df (10) 

Applying long wave length approximation {S ^ 1) and assuming the Reynolds number to small, 
the equation (8) approximates to 

i^yyyy - ot^'4^yy = (11) 

where 

.= =J^ + M (12) 
The solution of equation (11) satisfying (9) is 

^ = Ci + Csy + Cae"^ + C4e-"^ (13) 
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with 

^ _q^^ ahi{q + hi-h2) ^ (g + /ii - /i2)(e"('^i-''2) - 1)(1 - e-°^('^i-'*2)) 

^ 2 ^ a{hi — — 2tanh^{hi — A(l — cosha{hi — h2)) + 2a{hi — h2)sinha{hi — 

asinha{hi — /i2)(? + hi — /i2) 



C2 = -1 + 

C4 = 



2(1 — cosha{hi — /12)) + a(/ii — h2)sinha{hi — h2) 

{q + hi - /i2)e-"''i(e"('*i-''^) - 1) 
4(1 — cosha{hi — h2)) + 2a{hi — h2)sinha{hi — /i2) 
{q + hi- /i2)e"''2(e"(^i-^2) _ 1) 



4(1 — cosha{hi — /i2)) + 2a{hi — h2)sinha{hi — /12) 
where h2 < y < hi. 

In the absence of the magnetic field, these results are found to coincide with those obtained in 
[8], where the porosity effect is disregarded. 

The flux at any axial station in the flxed frame is found to be given by 



Q = / {u + l)dy = q + hi - h2, 

J /l2 



while the expression for the volumetric flow rate over one period T[— ^) of the peristaltic wave 
is obtained as 

Q = ^ f Qdt^^ [ (q + hi- h2)dt ^q + l + d 
I Jo I Jo 

The pressure gradient obtained from the dimensionless momentum equation for axial velocity 
can be expressed as 

rfp _ , _ 2 / _ 2 

Yyyy Yy 

a^sinha{hi — h2){q + /ii — /i2) 
2(1 — cosha{hi — h2)) + a{hi — h2)sinha{hi — h2) 
On integration over one wave length it reduces to 

Ap = -a^QIi - a^h 

where 



sinha{hi — h2)dx 



/o 2(1 — cosha{hi — /12)) + Oi{hi — h2)sinha{hi — h2) 
In the flxed frame of reference, the expression for the axial velocity is 



2(1 — cosha{hi — /12)) + 0({hi — h2)sinha{hi — h2) 
1 sinha{hi — h2){hi — h2 — I — d)dx 



U{X,Y,t) 



a{q + hi- h2){cosh^{hi - /12) - cosh^{hi + h2 - 2Y)) 
a{hi — h2)cosh^{hi — /i2) — 2sinh^{hi — /i2) 
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with 



hi — 1 + a cos[2n{X — t)] and /i2 



d - b cos[27r{X - t + (f))]. 



3 Results and Discussions 

Analytical expressions for the pressure gradient, the volumetric flow rate, axial velocity and 
the stream function for the problem under consideration have been calculated and presented in 
the previous section. In this section, the problem has been investigated further by employing 
appropriate computational method and numerical techniques. For the purpose of computation, 
the following set of values (which are in agreement with experimental data for physiological 
fluids) of the different parameters involved in the forgoing theoretical analysis have been used: 
a, b= 0.3, 0.4, 0.5, 0.7, 1.2; d=1.0, 2.0; Q = 1.4; M=0 to 8; /3k=0.02, 0.05, 0.08, 0.1, 0.13, 0.15, 
0.2, 0.5, 1, 10, 100, 1000; = 0.0, f , f , ^,7r. Small values of the permeability factor refer to 
cases when the medium is almost rigid. 

Due to the complexity of the problem under consideration, it was not possible to determine 
the analytical expression of Ap in terms of Q and other parameters. Thus had to resort to the 
use of appropriate numerical method along with quadrature formula. Similarly, an analytical 
treatment was found inadequate for finding the velocity as a function of Ap, X, Y, t and other 
parameters that are involved. This necessitated calculating first 



which was also a very complicated task. Here too, we had to make use of numerical quadrature. 
Velocity Change 

The distribution of axial velocity in the case of free pumping for different values of the phase 
difference, the permeability parameter k and the magnetic parameter M is shown in Figs 2-5. 
Fig. 2 shows that the results computed on the basis of our study in absence of any external 
magnetic field (M=0) perfectly match with the results reported in [8]. Figs. 3 (a)- (c) give the 
changes in the velocity distribution when M=2, as the phase difference (j) changes from to tt. 
These figures indicate that the velocity changes considerably with change in the values of the 



Ap + a^h 
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amplitudes of the peristaltic wave at the upper and lower walls of the channel. These figures 
reveal that the magnitude of the velocity gradually decreases and the extent of flow reversal near 
the boundary increases as the phase difference increases, when the wave amplitude at the upper 
boundary of the channel is more than or equal to the wave amplitude at the lower boundary. 

It is interesting to note that when the wave amplitude at the lower boundary exceeds that at 
the upper boundary, the scenario is similar when changes from to ^. In this case, when 
4> exceeds |, the flow reversal at the boundary starts decreasing and reduces to zero when 
When (p increases further, the flow reversal takes place in the central region and the velocity 
near the boundary increases. 

Fig. 4 shows that flow reversal takes places for all the values of the magnetic parameter exam- 
ined here and also that while in the central portion of the channel, the axial velocity decreases 
with the increase of magnetic intensity, the behaviour is slightly opposite in the lower and upper 
regions. It may be noted from Fig. 5 that in central region of the channel the velocity increases 
with the increase in the permeability parameter of the channel, while an opposite trend is found 
in the lower and upper regions. The occurrence of reverse flow at the initial stage is found to 
be independent of the permeability factor. 

Pumping Characteristics 

Figs. 6-8 illustrate the variation of the volumetric flow rate of peristaltic wave with pressure 
gradient for different values of the phase difference, permeability factor and the magnetic inten- 
sity. One may observe from Fig. 6 that in the range of values of the pressure gradient examined 
in the present study, in the entire pumping region (Ap >0), in the free pumping zone (Ap =0) 
as well as in the co-pumping region (Ap <0) for Ap >-12.3, the volumetric flow rate decreases 
with the increase in the phase difference. However, the trend reverses as soon as the pressure 
gradient drops below -12.3 . These observations, of course, hold good for a particular set of val- 
ues indicated in Fig. 6 caption. Fig. 7 shows that in the entire pumping region the volumetric 
flow rate decreases with the increase in permeability parameter, where as in the co-pumping 
region, a reverse trend is noticed. One can also observed that the volumetric flow rate can be 
gradually increased in the pumping region and gradually reduced in the co-pumping region by 
increasing the intensity of the applied magnetic field. 

Trapping 

A theoretical study on peristaltic pumping of a fluid was conducted by Shapiro et al. [31] for 



9 



the case when the wave length is large and the Reynolds number is small. This study revealed 
that at high flow rates and large occlusions, there is a region of closed stream lines in the wave 
frame, whereby some fluid gets trapped within a wave of propagation. The trapped fluid mass is 
found to move with a mean speed equal to that of the wave. Thus the phenomenon of trapping 
may be looked upon as the formation of an internally circulating bolus of fluid. Owing to the 
trapping phenomenon, there will exist stagnation points, where both the velocity components 
of the fluid vanish in the wave frame. The stagnation points will be located at the points where 
the central line of the channel intersects the curve ifj = {ijj being the stream function). In case 
of asymmetric channel flow, however, some stagnation points may be located elsewhere. 
Figs. 9(a)-(d) reveals a very interesting fact that with the same amplitude, a=b=0.5 for Q — 1.4, 
the bolus appearing in the central region for = moves towards left and decreases in size as 
(j) increases. For (j) — n, the bolus completely disappears and the stream lines are parallel to one 
another. 

The influence of the intensity of the applied magnetic fleld on the streamlines for peristaltic 
transport velocity where a=b=0.5 and /3k=1.0, (j)=0.0 and d=1.0 is demonstrated in Figs. 10(a)- 
(d). These figures indicate that when the magnetic fleld intensity is low, bolus is formed at the 

central region of the channel. These bolus is shifted towards the upper boundary of the channel 
with the increase in the magnetic field strength. Fig. 10(d) reveals that the bolus completely 
vanishes when M >7.0. 

Figs. ll(a)-(d) indicate that when M=l, in the absence of phase difference, bolus is formed for 
greater values of the permeability factor but no bolus is formed when the permeability parameter 
of the channel is small. 

A comparison between Fig. 11(d) and Fig. 12(a) reveals that when M=l, bolus formation is 
of different type when the phase difference is |, even the channel has the same permeability. 
But when the intensity of the magnetic field is higher, the bolus disappears even if the phase 
difference and permeability remain unchanged [cf. Fig. 12(b)]. Figs. 13(a)-(b) give us an idea 
of the type of bolus formation for greater values of the permeability factor. The variation in 
the shape of the streamhnes and bolus formation for different sets of values of a and b can been 
seen from Figs. 14 (a)- (c). 
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4 Summary and Conclusion 



The peristaltic transport of a Newtonian fluid (as an idealization of a physiological fluid) in a 
porous asymmetric channel under the action of an externally applied magnetic fleld has been the 
concern of the present investigation. The study particularly pertains to the situation when the 
Reynolds number is low and curvature of the channel is quite small. Emphasis has been paid to 
study some important bio-fluid dynamical phenomena of peristaltic transport such as the pump- 
ing characteristics, axial velocity as well as streamline patterns and trapping as a function of 
asymmetric wave motility parameters, the permeability parameter and the magnetic parameter. 
All these have been treated both analytically and numerically. The results and observations are 
good agreement with those reported in [8, 9, 10, 32]. 

On the basis of the observations of the present study, one can make an important conclusion 
that it is possible to increase pumping for the adverse pressure gradient by applying an external 
magnetic field and that the bolus formation can be avoided by suitably adjusting the magnetic 
field intensity. The present study reveals further that the velocity change is strongly dependent 
on the phase difference, amplitudes of wave propagating in the upper and lower regions, mag- 
netic field intensity and the channel permeability. 

The results throw some light on problems associated with fiuid movement in the gastrointestinal 
tract, intra- uterine fiuid motion induced by uterine contraction as well as fiow through small 
blood vessels and intra-pleural membranes. 

The study suggest that a convenient criterion for the presence of a trapped bolus is the existence 
of stagnation points in the wave frame. This physical phenomenon can be considered responsible 
for formation of thrombus in blood and the movement of food bolus in the gastrointestinal tract. 
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Fig. 1: A physical sketch of the problem 




Y ► 

Fig. 2: Distribution of axial velocity for X=l, t=l, a=0.7, b=0.7, d=l, (j)=0, /3k=1000, Ap=0, 

M=0 
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Fig. 3a: Distribution of axial velocity for different values of 0, where X=l, t=l, Ap=0, a=0.5, 

b=0.3, d=l, /3k=0.8, M=2 




Fig. 3b: Distribution of axial velocity for different values of 0, where X=l, t=l, Ap=0, a=0.5, 

b=0.5, d=l, /3k=0.8, M=2 
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Fig. 3c: Distribution of axial velocity for different values of 0, where X=l, t=l, Ap=0, a=0.3, 

b=0.5, d=l, /?k=0.8, M=2 
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Fig. 4: Distribution of axial velocity for different values of M, where X=l, t=l, Ap=0, a=0.5, 

b=0.3, d=l, ;3k=0.8, 0=f 
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Y ^ 

Fig. 5: Distribution of axial velocity for different values of /3k, where X=l, t=l, Ap=0, a=0.5, 

b=0.3, d=l, M=2, (/)=f 




Fig. 6: Variation of Q with for different values of <p along with a=0.7, b=1.2, d=2, /3k=0.1, 

M=2 
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Fig. 8: Variation of Q with Ap for different values of M along with a=0.7, b=1.2, d=2, 
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(d) 

Figs. 10: Streamlines when a=b=0.5, = 1.0, (f) = 0.0, d=1.0 (a) M=0.0, (b) M=3.0, 

(c) M=5.0, (d) M=7.0 
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Figs. 14: Stream lines when M=2, /3k=0.2, 0=f , d=l (a) a=0.3, b=0.5 (b) a=0.4, b=0.4 

(c) a=0.5, b=0.3 
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